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Abstract 

We study the evolutionary dynamics of an asexual population of nonmutators and mutators on a 
class of epistatic fitness landscapes. We consider the situation in which all mutations are deleterious 
and mutators are produced from nonmutators continually at a constant rate. We find that in an 
infinitely large population, a minimum nonmutator-to-mutator conversion rate is required to fix the 
mutators but an arbitrarily small conversion rate results in the fixation of mutators in a finite 
population. We calculate analytical expressions for the mutator fraction at mutation-selection 
balance and fixation time for mutators in a finite population when the difference between the 
mutation rate for mutator and nonmutator is smaller (regime I) and larger (regime II) than the 
selection coefficient. Our main result is that in regime I, the mutator fraction and the fixation time 
are independent of epistasis but in regime II, mutators are rarer and take longer to fix when the 
decrease in fitness with the number of deleterious mutations occurs at an accelerating rate 
(synergistic epistasis) than at a diminishing rate (antagonistic epistasis). Our analytical results are 
compared with numerics and their implications are discussed. 

KEY WORDS: mutators, genetic drift, epistasis, fixation time 
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Mutators are those ce lls or genotypes that have a higher mutation rate than the wild type (jVisser 



2002; 



Baeretal 



20071 ). While transient mutators arise and decline in response to an external 



perturbation such as stress, constit utive mutators carry defect (s) in p roofreading and mismatch 



repair pathways and are heritable (jMiller 



1996 



Rosenberg et al 



199a ). The latter class of mutato rs 



which we consider here 



rave been seen in several microbial 



1996; 



Matic et al 



experimental (jMao et al 



1997 



Oliver et al 



1997 



2000 



Biorkholm et al 



Sniegowski et al 



1997; 



jopulations in natural 



2001 



Boe et al 



Trobner and Piechocki 



fequency is seen to vary widely from less than a percent 



19971 ) 



1984; 



Mao et al 



Richardson et al 



LeClerc et al. 



2002 ) and 



200 0]) settings. The mutator 



Gross and Siege! 



19971 ) to hundred percent (jMao et al 



1981 



1997; 



Sniegowski et al. 



As most mutations are deleterious ( Drake et al. 



19981 ). how do mutators manage to reach significant 



frequencies? In a maladapted population, while the mutators produce more deleterious mutations 
per generation compared to the nonmutators, they also give rise to rare beneficial mutations more 



often. In an asexual popu 
population can hitch hike ( 



eventually get fixed 



many theoretical 



ation in which the mutator and its effects remain linked, the mutator 



Smith and 



Sniegowsk 



Taddei et al 



And 



19971 : 



et al 



1997; 



aigh 



1997 



19741) with benefici al mutations to high frequencies and 



Shaver et al 



Tenaillon et al 



e and Godcl ld 120061) I and experimental studies (ITrobner and Piechocki 



Cooper and Lenski 



2000 



Shaver et al. 



2002 



199 



2002 1). This situation 



las be en a subject of 



Pal mer and L ipsitch 



Gentile et al 



2011 



1984 



2006 



Sniegowski et al 



Here we are interested in the dynamics of an asexual population of nonmutators and mutators after 
the adaptation process is over. We assume that only deleterious mutations can occur and neglect 
rare beneficial mutations aris ing due to back or compensatory mutations and the creation of 



nonmutators from mutators ( Drake et al 



199a ). In this setting, clearly the mutators cannot go to 
fixation by hitchhiking. However if the mutators are being continually produced and the population 
size is finite, it is obvious that they will eventually get fixed for arbitrarily small conversion rate. It is 
important to note that in an infinitely large population at mutation-selection balance, a minimum 
nonmutator-to-mutator conversion rate is required above w hich the mutators get fixed 



(Tannenbaum and Shakhnovich 



Some recent works ( 



2005; 



Nagar and Jain 



Soderberg and Berg 



2011 



Lynch 



2009). 



2011 



have studied the dynamics of mutator 



fixation in finite populations when beneficial mutations are absent assuming that the fitnesses do not 
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interact epistatically and the mutational effects are w eaker than the selective co s ts. However 



experiments i ndicate that epistasis is rather common (jde Visser and Elena 



'hillips 



Jain 



2008, 



2007 



Kouvos et al 



2007; 



20081) an d theoretical studies have shown that it strongly affects both statics and dynamics 



20101 ) . To our knowledge, the effect of epistatic interactions on mutator dynamics has 
not been explored in previous theoretical studies. On the experimental front, strong epistatic 
interactions betwe en deleterious mutations ha ve been observed in S. typhimurium for a broad range 



of mutation rates ( Maisnier-Patin et al. 



2005|). 

In this article, we consider the evolution of both infinitely large and finite sized populations on a 
class of epistatic fitness landscapes. The epistatic interaction between the fitnesses is antagonistic 
(synergistic) if the fitness of the double mutant is larger (smaller) than the sum of the fitnesses of 
the two single mutants. Using a continuous time, deterministic model, we first show that in an 
infinite population, a transition occurs at a critical conversion rate between a mixed state with both 
nonmutators and mutators and a pure mutator state for any epistasis. We then develop a diffusion 
theory to calculate the time to fix mutators when the conversion rate is small. We find that the 
fixation time depends weakly on the population size for small populations but grows exponentially 
fast for larger populations. The population size at which this crossover occurs and the growth rate of 
fixation time for large populations are governed by the mutator fraction at mutation-selection 
balance which we analytically calculate when mutational effects are weaker (regime I) and stronger 
(regime II) than the selective effects. We find that in regime I, the mutator fraction is independent 
of the epistasis parameter (see (I24p ) and as a consequence, the fixation time is also insensitive to 
epistasis. But in regime II, the mutator fraction exhibits a nontrivial dependence on epistasis (see 
(I27|) ) and decreases as epistasis changes from antagonistic to synergistic. As a result, the fixation 
time for synergistically interacting fitnesses is larger than that for antagonistic ones. The 



implications of our results for the evolution of mutation rates are discussed. 



Models 

We consider a haploid, asexual population evolving on a fitness landsca pe in w 



rich the Malthusian 



fitness F(k) of a genotype carrying k deleterious mutations is given by (IWiehe 



19971 ) 



F(k) = -sk a 
4 



(1) 



where s is the cost of deleterious mutation and the parameter a > controls epistasis. If a = 1, 
there is no epistasis while a < 1 and > 1 respectively correspond to antagonistic and synergistic 
epistasis. A deleterious mutation occurs at a rate U for a nonmutator and V = XU, A> 1 for a 
mutator. Besides replication and mutation, a mutator is created from a nonmutator at a forward 



convers ion rate f 



b <C / (jNinio 



j ut the reverse re action is ignored since such events occur at a backward rate 



1991 



Boeetal 



12000). 

We implement the stochastic dynamics of a population of fixed size iV using the standard 
Wright-Fisher process in which each individual in the new generation chooses a parent with a 
probability proportional to the parent's fitness. Then the mutations chosen from a Poisson 
distribution with mean given by the parent's mutation rate are introduced and finally a nonmutator 
offspring is converted into a mutator. 

In an infinite population in which the time evolution occurs deterministically, the average fraction 
V{k, t) and Q(k, t) of nonmutators and mutators respectively with k deleterious mutations at time t 
evolves according to the following coupled equations: 



dV(k,t) 

di 
dQ(k,t) 

dt 



U(V(k - 1, t) - V(k, t)) + (F(k) - W(t))V(k, t) - fV(k, t) 
V(Q(k - l,t) - Q(k,t)) + (F(k) - W(t))Q(k,t) + fV(k,t) 



(2) 
(3) 



In the above equations, the fractions V(—l,t) = Q(—l,t) = and 

oo 

W(t) = Y J F{k){V{k,t) + Q{k,t)) 



(4) 



fc=0 



is the population fitness which ensures that the normalisation condition Yl^o^i^^) + Q(k, t) = 1 is 
satisfied in all generation s. The above equations can be derive d using standard meth ods 



(jBaake and Gabriel 



2000) starting from a discrete time model (jNagar and Jain 



20091 ). 



Results 

FIXATION OF MUTATORS IN INFINITE AND FINITE POPULATIONS 

We first show that in an infinite popul ation at mutation-selection balance, a minimum conversion 



rate f c is required to fix the mutators (jTannenbaum and Shakhnovich 



2005: 



Nagar and Jain 



20091 ) 



Although the nonmutator population carries less mutational load than the mutators, it suffers a loss 



due to the creation of mutators when the rate / is nonzero. As a result, there is a transition at a 
critical conversion rate f c below which the population is in a mixed state composed of both 
nonmutator and mutator subpopulations and above which only mutators are present (see inset of 

Fig. nj. 

To see this, consider the nonmutator and mutator frequencies at mutation-selection balance when 
they become time-independent. Then ([2]) and ([3]) reduce to 



UV(k - 1) + (F(k) - W - f - U)V(k) = 
VQ(k-l) + (F(k)-W-V)Q(k) + fV(k) = 



(5) 
(6) 



Mixed state: It is clear from ([5]) that a nonzero distribution of nonmutators is possible only if the 
fraction 'P(O) is nonzero. For k = in (J5]), we immediately have 



-f-U , 7>(0) ^0 



(7) 



Thus the mean fitness in the mixed state is independent of the epistasis exponent a, mutator rate V 



and selection effect s. T 



lis result is a simp le generalisation of the well known result for / = 



( Kimura and Maruvama 



1966; 



Haigh 



19781 ) 



Pure mutator state: If the nonmutator frequency in the fittest class is zero, ([6]) reduces to 



VQ(k - 1) + (F(k) -W- V)Q(k) = 



(8) 



where a tilde is used to refer t o quantities in t 



have been studied extensively (jJain and Krue 



shown that an error threshold transition ( 



re pure mutator state. Equation (E l) and its relatives 



2007: 



Eigen 



Jain and Seetharaman 



2011 



) and it has been 



197ll ) can occur in which a population localised 



around the fittest sequence gets delocalised beyond a critical mutation rate. For the fitness function 
(PP), such a transition occurs in the special case of sharp peak fitness land scape (a — > ) but the 



Baake and Gabriel 



population remains in th e localised state at any mutation rate for a > (jWiehe 
"J. 



19971 : 



2000). Since the population fraction in the fittest class is nonzero in the localised 



state, using (E| for k = 0, we get 



-V , Q(0) / 



(9) 



Critical conversion rate: As shown in the inset of Fig. [lj with increasing rate f , the mixed sta te 



fitness W decreases linearly until it reaches the pure mutator state fitness W (jNagar and Jain 



2009|). 



Matching the mean fitnesses (J7|) and ([9]) in the two states gives the critical forward conversion rate as 

fc = V-U (10) 

which does not depend on the epistasis parameter a. (However for sharp peak fitness landscape, an 

additional critical conversion rate exists besides (|10j) . which is discussed in Appendix A.) 

While a minimum conversion rate is required to fix mutators in an infinitely large population, the 

mutators get fixed for any / > when the population is of finite size. Starting with a population of 

N nonmutators each with fitness one, we follow the dynamics of the total nonmutator frequency 

P(t) = J]fc=o P(k, t) where P(k, t) is the instantaneous fraction in the fitness class k. As shown in 

Fig. [TJ the fraction P(t) first relaxes to a value close to the deterministic nonmutator fraction 

*P = /Cfclo^^) (phase 1) about which it fluctuates (phase 2) until it enters the absorbing state with 

zero nonmutators (final state). 



TIME TO FIX MUTATORS IN A FINITE POPULATION 

Ignoring the time for the nonmutator population to relax to its deterministic frequency, we compute 
the time for the loss of nonmutators starting from deterministic fraction V using a diffusion theory 



( Ewens 



19791 ). We begin by calculating the mutator fraction at mutation-selection balance as it is 



required to find the fixation time. In the following, we work in the biologically relevant limit / — > 



( Ninio 



1991 



Boe et al. 



2000j) and also assume that the epistasis exponent is nonzero and finite 



(0 < a < oo). 



Mutator fraction at mutation-selection balance 

To calculate the mutator fraction, we first expand the population fractions in a power series about 
/ = by writing 

fn Vn{k) 



r(k-j) = Y,f r 



71=0 

oo 



Q(k;f) = £/ 



nQn(k) 



n=0 



)V. 



(11) 

(12) 



where P n {k) and Q n (k) respectively denote the nth derivative of the nonmutator and mutator 
fraction with respect to rate / evaluated at / = 0. Since Qo(^) = Q(k; / = 0) = for all k 



7 



(jNagar and Jainll2009l ). the mutator fraction Q(k; f) ~ fQi(k) where we have neglected terms of 
order f 2 and higher in (fT2|) . On using ([7J), (fTT]) and (fT2|) in the steady state equation ([6]), a 
straightforward calculation shows that Qi(k) obeys the following equation: 

VQi(k - 1) + (F(k) - 5U)Qi{k)+V {k) = (13) 

where 5U = V — U. Similarly using j5]),([7|) and (jlip . the equation for nonmutator fraction Vo(k) can 
be written as 

UV (k - 1) + F(k)V (k) = (14) 
which can be readily iterated to yield 



V (k) = - 



uyv (o) 



(kiy 



(15) 



where the fraction Vq(0) is given by 



V (0) 



n -1 



3=0 



(16) 



The sum on the right hand side (RHS) of the above equation does not appear to be kno wn in term s 



of standard mathematical functions except for some special values of a such as 1 and 2 (|Jain 
For the mutator fractions, (I13D gives 

O (k) = V ^ k ~ l ) + ^o( fc ) 
^ U ' SU-F(k) SU-F(k) 

On using Qi(— 1) = in the above equation for k = 0, we get 

Po(0) 



20081 ) 



Si(0) 



<5£7 



(17) 



(18) 



For other fitness classes, on repeatedly iterating (fTT|) . we find that 



Qi(fe) = Qi(o) ] 



v 



r 



. 5t/ - F(m) ' 



A; > 



(19) 



J- c)£7 — F(m) V ^ 

Finally on summing over all the fitness classes, we obtain the mutator fraction to leading order in / 

as 



oo k 



fe=0 j=0 



g(fc) / 

G(i) 5C/ + sj< 



(20) 
(21) 



where we have defined 



G(k) = n 



V 



m=l 



SU + sm a 



(22) 



To evaluate the mutator fraction in (|2ip . we first estimate the product G(k) as explained in 
Appendix B. Using the result (|5ip in (|2ip . we find that 



Q~su^^{5u) exp 

k=0j=0 



s (k a+1 -j a+1 ] 



VoU) 



(23) 



SU a + 1 

Regime I: If s/SU 3> 1, the dominant contribution to the inner sum on the RHS of (|23p comes from 
the term with j = k. Neglecting the summand with j ^ k, we immediately obtain 



Q « , «5[// s « i 



(24) 



Note that this solution is independent of the selective cost s and the epistasis exponent a. 
Regime II: For s/SU <C 1, the evaluation of the double sum in (|23p is nontrivial and we calculate the 
mutator fraction in this regime within a Gaussian approximation as explained in Appendix C. 
Plugging (|52|) and ([56]) in the expression (I2ip for mutator fraction and approximating the sums by 
integrals, we get 

aU{k-k*) 2 f k aUU-k*) 2 a(j-k*) 2 



f 



Q ~ T77 



517 V 2vrfc* jo 
/ 



e" 



erf(V/c*)+erf ( x j — (k - k* 



a 



25U jo 

where erf (a;) is the error function. For A;*, A S> 1, the second integral on the RHS of the above 
equation vanishes and the first integral can be readily evaluated to give 



(25) 
(26) 



Q 



For a = 1, we recover the known results (jNagar and Jain 



7T U- 



\ 2a (A-l)s^ 



dU/s » 1 



(27) 



2009 



Desai and Fisher 



2011 



Our analytical results (J24j) and (|27p are compared with the exact expression (I2ip evaluated 
numerically in Figs. [2][1] and show a good agreement. A detailed discussion of these results is given in 
the last section of this article. 



Diffusion approximation for fixation time 



We now develop a diffusion theory (jEwens 



19791 ) to estimate the average time to lose all the 



nonmutators. Using a backward Fokker-Planck equation, it can be shown that the average 



absorptio n time T to t 



given by (|Ewens 



1979; 



le state w ith zero nonmutators starting from total nonmutator fraction V is 



Jain 



2008) 



dx 



2 / dy i{j(y) [ 

JO Jy D 2 {x)llj{x) 



(28) 



where the function 

^( x ) = e J (29) 
and Di(x), D2(x) are drift and diffusion coefficients respectively. 

To calculate an expression for the drift coefficient D\, we first sum over all k > in (|2|) and find that 

&P(t) 



Of 



W P (t)-(f + W(t))V(t) 



(30) 



where V(t) = YlT=o'^ > (^^) an d VVp(f) = Sfclo F{k)'T ) {k, t). At mutation-selection balance, due to 
([7]) and ([30]) . the mean fitness of the nonmutator population is given by 



W P = -UV 



(31) 



Using the definition of drift coefficient ([Ewens 



1979 



Jain 



20081 ) and the result ([30)) . we may write 



D 1 (P(t)=V(t)) 



dV(t) 

at 

D ( WP t W 

Px l—-f-w 



(32) 
(33) 



Since the above eq uation does not close in P(t), we use a 



expression for D\ (jStephan et al. 



1993 



Stephan and Kim 



simp le argument to find an approximate 



20021 ) . We first note that D\ has two zeros 



namely and V since the nonmutator fraction P{t) does not change with time in the absorbing state 
and (quasi)steady state. Therefore we can write 



D, 



(34) 



where C is a proportionality constant which can be determined by considering the second factor on 
the RHS of (I33p when the average population fraction is close to the mutation-selection balance 
(phase 2). Assuming that the population fitnesses Wp{t) and W(t) have the same functional form in 
phase 2 as in the time-independent steady state, we can express them in terms of P(t) as 



W P (P) 
W(P) 



-UP 



-u 



1-p 



(35) 
(36) 



10 



where we have used (|3ip and eliminated / from ([7]) in favor of V using Q = fQ\. A power series 
expansion of the preceding equations in P — V gives 



W P 
~P~ 



f-W 



Qi 



f 



1 



P 
V 



(37) 



On using the last expression in f|33[) and comparing it with (|34f) . the constant C is immediately 
determined. 

Furthermore as the number of offspring produced per generation are distrib uted according to a 



binomial distribution, the variance in the offspring number gives (jjain 

P(l - P) 



2008) 



Do 



N 



(38) 



Using the above expressions for D\ and D2 in (|28p . the fixation time for mutators can be estimated 
as described in Appendix D. The final result is 



2WJ 1 



2Nf 

e 2 — e 



27V/ ) , JV « 1/(2/) 



, 2 2NfV 

>e— , N » 1/(2/) 



(39) 



2Nf A V 2 

where the mutator fraction Q is given by ()24p or (|27p . The fixation time obtained using the above 
expressions and numerical simulations is shown in Figs. [5] and [6] for representative values of a. The 
following section discusses these results in more detail. 



Discussion 

In this article, we studied the dynamics of a population of nonmutators and mutators, starting from 
a well adapted population of nonmutators, when the mutations are unconditionally deleterious. 



EFFECT OF GENETIC DRIFT 

The nonmutator population decreases due to a conversion into mutator genotypes. But in an 
infinitely large population, this loss can be compensated if the conversion rate is sufficiently small. 
On a class of fitness landscapes defined by ([1]), a transition occurs at a critical rate f c between a 
mixed state with nonmutators and mutators and a pure mu tator state 



( Tannenbaum and Shakhnovich 



2005; 



Nagar and Jain 



2009). Interestingly, the critical conversion 



11 



rate given by (fTUl) is indepen dent of a (but see Appendix A also). As the previous studies (jjohnson 



1999; 



Desai and Fisher 



2011 



) on infinite populations worked exc 



usively in the sma 


11 f limit, the 


Drake et al. 


1998; 


Boe et al. 


2000) 



and A = 10 in (|10p . we find that the transition occcurs at f c ~ 10 2 which is several orders of 
magnitude larger tha n the spontaneous conversion rate / ~ 10 -6 — 10 -7 estimated in E. coli (]Ninio 



1991 



Boe et al 



2000) 



The difference between an infinite and finite population is that in the latter case, the nonmutators 
get wiped out for arbitrarily small / in a finite time as illustrated in Fig. (TJ If the con version rate i s 



set to zero in a finite population of nonmutators and mutators, either of them can fix (|Wvlie et al 



20091 ) . But for nonzero /, due to the simple fact that the mutators are being continually generated, a 



finite population of nonmutators will be converted into mutators in a finite time. 



EVOLUTIONARY REGIMES 

The mutator fraction at mutation-selec tion balance ha s been calculated by seve ral authors on 
non-epistatic fitness landscapes. While 



Johnson 



(1999); 



Desai and Fisherl (120111 ) computed it for a 



general distrib ution of fitness e f fects and small conversion rate /, an exact solution was obtained for 



arbitr ary f by 



work, 



Nagar and Jain 



Desai and Fisher 



2011 



20091 ) when all mutations have the same selective effect. In their 
have pointed out that the behavior of the mutator fraction depends 
on the relative strength of mutational effects (SU) to selective effects (s). When selective 
disadvantage exceeds the mutational e ffects, the mutations can be treated as effectively lethal 



( Johnson 



1999 



Desai and Fisher 



20111 ). This is because when SU <C s (regime I), the most 
populated nonmutator fitness class k* in (|48j) is below unity and therefore the nonmutator 
population is effectively localised in the fitness class with zero deleterious mutations. As a 
consequence, mutators which arise due to nonmutators, are also localised in the zero mutation class. 
Setting all the nonmutator and mutator fractions other than in the zeroth mutational class to zero, 
(118|) immediately leads to (|24p which is (unsurprisingly) independent of epistasis. Moreover the 



rcsu 



2011 



t (1241) is expected to be insensitive to the distribution of selective effects also (jDesai and Fisher 
) since both subpopulations are localised in the zero mutation class. However in the opposite 
limit SU 3> s (regime II), mutator subpopulation is spread over many fitness classes and the simple 
argument given above underestimates the mutator fraction as can be seen by extrapolating the 



12 



solution (|24p to regime II in Figs. [2 M1 An expression for mutator frequency in regime II has bee n 



obtained in the absence of epistasis 



Johnson 



1999; 



Nagar and Jain 



2009 



Desai and Fisher 



2011 



and here we have extended the previous analyses by including it. 

Figures [2]|3] show that the mutator fraction in regime II decreases monotonically as the epistasis 
parameter a increases. To understand this behavior, we first note that the fitness difference 
F(k) — F(k + 1) ~ sak a ~ l for large k approaches zero for a < 1 and infinity for a > 1 as k increases. 
Thus the mutators stand a better chance of survival when the fitness interactions are antagonistic 
than synergistic. For a given a, as large selective cost is detrimental for the mutators, the mutator 
frequency decreases (regime II) to a constant (regime I) as s increases (refer Fig. [2]) . The effect of 
increasing mutator strength A is also to reduce the mutator fraction as shown in Fig. [U The 
variation of mutator frequency with the spontaneous mutation rate U shown in Fig. [3] however 
displays an interesting feature: as U increases, the mutator frequency decreases monotonically for 
a > 1/2 but for a < 1/2, it decreases in regime I and increases in regime II. At a = 1/2, the mutator 
fraction decreases to a constant as mutation rate U is raised. This behavior may be understood as 
follows: in regime II in which both subpopulations are spread over the fitness landscape and carry 
many deleterious mutations, both nonmutators and mutators incur a negligible fitness cost when a 
mutation occurs for small a but due to the unidirectional flux from nonmutators to mutators, the 
mutator fraction approaches unity as U increases (also see Appendix A). This suggests that even if 
the spontaneous mutation rate of the nonmutator is high, the mutators can still thrive provided the 
fitness cost increases sufficiently slowly. 



On non-epistatic fitness landscapes, the short time dynamics ([Johnson 



1999; 



Desai and Fisher 



2011 



until the population reach e s a mutation -selection balance and the long term dynamics 



(jSoderberg and Berg 



2011 



Lynch 



2011 



in which the mutators get fixed due to stochastic 
fluctuations have been studied. Here we have developed a diffusion theory to find the time to fix 
mutators in a population of size N when the conversion rate / is small in the presence of epistatic 
interactions. Our results given in (|39[) predict that the fixation time remains roughly constant for 
small populations and increase exponentially fast beyond a critical population size N c ~ Q/ _1 . Thus 
in the limit N — > oo, the time to fix mutators is also infinite which is consistent with the result that 
an infinitely large population with conversion rate / -C f c exists in a mixed state. The fixation time 



depends on the biological parameters U,\,s,a through the mutator fraction at mutation-selection 
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balance. In regime I, since the mutator fraction is independent of a, this has the direct consequence 
that the fixation time is also unaffected by a (refer Fig. [5]) and therefore experiments aiming to 
detect epistasis using mutators should not operate in regime I. From (|24p . we find that the crossover 
population size iV c = (2JC/)" 1 is approximately 25 for the parameters in Fig. [5] in agreement with the 
simulation results shown in the inset. The fixation time in regime II depends on the epistasis 
exponent in a nontrivial manner. As discussed above, the mutator fraction in regime II decreases 
with increasing a and therefore due to (|39p . it is expected that the time to fix mutators increases 
and the crossover population size decreases as a increases. These expectations are borne out by the 
simulation results in Fig. [6j We also note that since the mutator fraction in regime II is higher than 
in regime I, the fixation time in regime II is shorter than that in regime I. Although the diffusion 
theory developed here using simple arguments captures the essential features of the fixation time, a 
more sophisticated analysis based on a path-integr al formulation may be requ ired for a better 



quantitative agreement with the numerical results (jNeher and Shraiman 



20121 ) 



EVOLUTION OF MUTATION RATES 

In this article, we studied the dynamics of the evolutionary process by which the mutation rate of a 
finite asexual population increases from U to XU in time T\ = T. This population can in turn 
produce mutators with mutation rate X2V at a rate fa which take an average time T2 to fix and so 
on. Thu s the mutatio n rate of the population can keep increasing in a manner similar to Muller's 



ratchet ( Muller 



1964! ) in which an asexual populatio n keeps accumulating delete rious mutat i ons. 



Soderberg and Berd (|201ll ) and 



Lynch! |201l|) 



This ratcheting process has been studied recently by 
when epistatic interactions are absent and mutator strength is weak (regime I). Here we have 
focused on the dynamics of the first click of the ratchet for a more general set of parameters. 
When epistasis is absent, we can apply the results obtained here for the first click of the ratchet to 
the succeeding clicks. Assuming that the conversion rate and mutator strength is same for all clicks 
of the ratchet, the mutation rate at the nth click of the ratchet will increase to U n = X n U. Then 
using (|24p and (|27p in the expression (J39]) for fixation time, we find that the time elapsed between 



(n — l)th and nth click grows as e NSU " if 5U n -i < s/X and e N ^ SUnS otherwise, where 

SU n = U n +i — U n . Thus the time between consecutive fixations increases but the rate of increase is 

slower in regime II than in regime I. Our analytical estimates above are consistent with the 
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numerical results of 



Soderberg and Berg] (|2011 



who observed that the fixation time increases with 



mutatio n rate U„, and mutator strength A. Our results (|24p and (|39p are also in agreement with 



those of 



Lynch 



2011 



who, using a diffusion theory and simulations, showed that the mean time to 
fixation is roughly constant in SU n and selective effects during the initial clicks characterised by 
SU n <C 1/N but increases exponentially fast for later clicks. On epistatic fitness landscapes, as the 
time between su ccessive clicks of the ratchet depends on the fittest class available at the end of a 



click (Jain 



20081 ). a quantitative understanding of the successive fixation process requires more work 
and will be discussed elsewhere. Here we merely speculate that for a > 1/2, every successive fixation 
event will take longer to occur since the mutator fraction decreases with increasing mutation rate 
and may result in fitness reduction to an extent that the population may not remain viable. On the 
other hand, as the mutator fraction at mutation-selection balance increases with mutation rate for 
a < 1/2, the ratchet is expected to click fast but the population may be able to tolerate the errors as 
the fitness penalty does not increase commensurately. Interestingly in an experiment on bacterium 
S. typhimurium measuring the decrease in fitness as deleterious mutations accumulate, the mutation 
rates was varied up to hundred times more than the wild type and the fitness function of the form 
( J) was found to be consist ent with the epistasis exponent awl/2 for different mutation rates 



(|Maisnier-Patin et al 



20051 ). Further experimental and theoretical studies of the mutator dynamics 



on epistatic fitness landscapes are required to better understand the ratcheting process. 



EFFECT OF BENEFICIAL MUTATIONS 

In the discussion so far, we have assumed that the backward conversion rate b at which a mutator is 
created from a nonmutat or is zero. T his assumption is motivated by experiments on E. coli which 



indicate that b/f ~ 10 3 ([Ninio 



1991 



j). Here we briefly discuss how a nonzero backward rate b affects 
our main results. An immediate consequence of the nonmutator creation is that the population 
remains in a mixed state, irrespective of its size. In the inset of Fig. the nonmutator fraction as a 
function of time is shown when the backward rate b is nonzero. We find that the nonmutator 
frequency is high for about 4000 generations after which the mutators take over. However as the 
nonmutators can be created at a rate b <C /, the mutators cease to dominate after 6000 generations 
and the process repeats. In Figs. [5] and [6l the results of our numerical simulations measuring the 
average time at which the nonmutator fraction becomes zero for the first time is compared with the 
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fixation time for mutators when 6 = for a wide range of parameters and we find that the loss time 
for nonmutators is mildly affected when the mutators are allowed to convert to nonmutators. 
However the time during which mutators dominate (data not shown) can not be obtained by merely 
interchanging the nonmutator and mutator labels since the mutators by definition carry a higher 
mutation rate. A complete analysis of this model is nontrivial and we hope to address it in a future 
study. 

ieie we focused on the mutator fixation when the mutations are deleterious but theoretical studies 



( Andre and Godelle 



2006 



Gerrish et al. 



2007; 



Wvlie et al 



2009 ) indicate that the mutators will take 



Ravnes et al 



(|2012j ). the number of 



over during adaptation also. In this scenario, besides the continual production of mutators and finite 
population size, an additional factor namely availability of beneficial mutations is at work. As a 
result, the mutators are likely to get fixed in a time shorter than the one computed here. Indeed our 
numerical results shown in Fig. [7] for a model in which bene ficial mutations can also occur are in 
agreement with this expectation. In a recent experiment by 
generations for the mutators to take over an adapting population of S. cerevisiae was measured. The 
fixation time was observed to increase with population size but in very large populations, the 
mutators were not seen to fix. The latter result is surprising since we have shown here that a finite 
asexual population will fix mutators even in the absence of beneficial mutations. As the fixation time 
is large for larger populations, a longer run of the said experiment may shed some light on the 
dynamics of mutator fixation during adaptation. 



Appendix A: Infinite population on sharp peak fitness landscape 

When the exponent a — > in ([I]), we get a sharp peak fitness landscape defined as 

F(k) = -s(l - 6 k , ) (40) 

The conversion rate at which the transition occurs depends on the nature of the pure mutator state. 
If the mutator population is in the localised state, (jlOp holds. But if the mutator population is 
uniformly distributed with an average population fitness W* = — s, we have 

ft = s - U , U < s (41) 
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Moreover, using (|40p in ([5]) and ([6]), we obtain 



V(k) = V(0) - , U < s (42) 



s 

l/-(7 + s-/ K - U + s(l - d fc ,o) - / 

On summing over k from zero to infinity on both sides of the above equations and using V + Q = 1, 
we get 

P= H) ('"D (44) 

which shows that the total nonmutator fraction vanishes at both f c and /* in accordance with the 
above discussion and the pure mutator state is the only solution for / > f c , f*. 



Appendix B: Evaluation of product (22) 



The product (|22p can be rewritten as 

G{k) : 



(m) exp 



v_ 

su 

V 

W 



exp 



exp 



sp 



S>( 1 + & 

o * ln ( 1+ ^ 

s 



dz ln(l + z c 



(45) 
(46) 
(47) 



where we have approximated the sum by an integral. The function G(k) is nonmonotonic in k and 
has a maximum at k = k* given by 

(48) 

This can be seen by considering the ratio G(k + 1)/G(k) which crosses unity at k = k* . As the 
population fraction T , o(k) given by (|15p also has a maximum at k*, the dominant contribution to the 
sums in (|21h comes when the argument k ~ 0(/c*). As a result, the upper limit of the integral on the 
RHS of (|4"T|) scales as 

i 

U 



5U 



< 1 



(49) 
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if the mutator strength A S> 1. 



This allows us to write 



G(k) 



V 

su 

v_ 

su 



exp 



exp 



SU^~ 

s k a+1 
SU a + 1 



KmY' 



dz z c 



(50) 
(51) 



Appendix C: Gaussian approximation in regime II 

To find the mutator fraction in regime II, we first note that both Vo^k) and G(k) have a maximum 
at k* given by (148 p . Using the Stirling's formula k\ ~ V2irk(k/e) k in p5j) . we approximate the 
fraction Vo(k) by a Gaussian centred about k* and obtain 



V (k) 



e 2k* 



2-rrk 



(52) 



which is normalised to unity for k* S> 1. The product (|5ip can also be approximated in a similar 
fashion and we have 

a k a+1 — k* a+1 



V 



(k-k*)\n( — 
(fc-F)ln 



5U J 5U a + 1 
V \ s k* a+1 



5U J 5U a + 1 
(k-k*) 



k-k*\ a+1 
l + ^—l 



k* 



aU 



*\2 



2k*5U 



{k-k*) 



(53) 
(54) 
(55) 



For strong mutators, we thus obtain 



aUtk-k*) 2 

G{k) w e 2E^i^G(r) 



(56) 



Note that the width of the above Gaussian approximating the product G{k) is about A times larger 
than that of the Gaussian in (1521) . 



Appendix D: Calculation of fixation time 

Using the drift coefficient (|34|) and diffusion coefficient (j38|) in the expressions (j28|) and (|29p . we get 



T = 2N I dy i>{y) / — -— - (57) 



jy X(1-X)1P(X) 
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where 



ip{x) 



2Nx 

ST(1 - x)~ 2Nf 



and V = 1 — fQi- Thus we have 



V 



_-lNjj_ 



2N I dy e Oi (1 _ y)~ 2A ^ / dx 
'0 Jy 



2Nx 

e^T(l-x) 2N f 



x(l — x) 

We first evaluate the integrals in the above equation for 2Nf < 1 by writing T = I\ + 12 where 

2Ny nltT , 2iV:c 



(58) 



(59) 



N2JV/-1 



2iV / dy e s i (1 - y)~ 2Nf / dx e s i (1 - x) 

JO Jy 

r"P 2Ny /"l 2JVa: 

2iV / dy e~^T(l -y)~ 2Nf / dx e"3T(l - x) 2 ^" 1 

j ./j/ 



(60) 
(61) 

Since the factor (1 — x) 2N ^~ l in the inner integral on the RHS of (|60p diverges at unity for 2N f < 1, 



the dominant contribution to the integral I\ comes due to this factor: 



V 



2Ny 



2N I dy e d (1 - y) 
Jo 

2l < e W- e ^f 



-2Nf 



f dx e^(l-x) 2Nf - 1 
Jy 



2Nf 



(62) 
(63) 

(64) 

(65) 
(66) 

Since 2Nf < 1, the behavior of the fixation time for N <C Qi/2 is dominated by (|64[) , Thus the 
diffusion theory predicts that the fixation time for small populations is a constant in N. We remark 
that the iV-independent behavior of fixation time will no t be obtained if instead of (j38j) . Feller 



If N <C Qi/2, the above integral gives the fixation time to be independent of N i.e. 

h-j-Qi = j,N<^ Qi/2 
On the other hand, for small 2Nf and 2N/Qi, the integral I2 may be approximated as: 

I 2 « 2N [ dy I' dx x~ l 

Jo Jy 

m 2NV , N < Qi/2 



diffusion D2 = P/N is assumed (see, for e.g., 



Jainl (|2008h ). 



For large populations with 2N f,2N/Qi 3> 1, a rough estimate for the fixation time may be obtained 
using (|59p . Approximating (1 — x) 2N f ~ e~ 2N ^ x , we get 



T « 2iV / dy e" — e 27V/2/ / 



dx e Q i e 



-2Nfx 



^1 g — 



2iV7> 2 

Thus for large iV, the fixation time increases exponentially in iV. 



(67) 
(68) 
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Figure 1: Loss of nonmutators in a finite population (main) due to stochastic fluctuations and continual 
production of mutators at a rate / = 10 -3 and in an infinitely large population (inset) at a critical 
conversion rate f c . The horizontal line in the main figure is the nonmutator fraction at mutation- 
selection balance. This fraction is shown in the inset on the left y-axis as a function of /. The inset 
also shows the mean fitness W (thick line) for an infinitely large population on the right y-axis. The 
other parameters used in all plots are s = 10~ 2 , U = 2 x 1CT 2 , A = 5, a = 1. 
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Figure 2: Variation of mutator fraction at mutation-selection balance with selective cost s for a = 
l/2(o), l(o), 2(A). The other parameters are / = W~ 7 ,U = 1(T 4 ,A = 10 2 . The exact sum fl2T} is 
shown by points and the asymptotic expressions (|24p and (|27p by solid curves. The inset shows the 
data in the main figure when plotted as a function of selective effect s. 
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Figure 3: Variation of mutator fraction at mutation-selection balance with mutation rate U for a = 
1/4(V), l/2(o), l(o), 2(A). The other parameters are / = 1CT 7 , s = 1(T 2 , A = 10 2 . The exact sum |(2TJ) 
is shown by points and the asymptotic expressions (|24p and (|27p by solid curves. The inset shows the 
data in the main figure when plotted against the exponent a at fixed 5U/s = 990. 
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Figure 4: Variation of mutator fraction at mutation-selection balance with mutator strength A for 
a = l/2(o), l(o), 2(a). The selective cost s is 10~ 2 (main), 10 -6 (inset) while the other parameters 
are / = 10 _7 ,C7 = 10 -4 . The exact sum (|2ip is shown by points and the asymptotic expressions 
(1241) and (I27p by solid curves. The result (I27j) is not shown in the main figure as it holds when 
k* = (U /s) 1 /" s> 1 but it is compared with the exact expression (|2ip in the inset. 
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N 

Figure 5: Variation of fixation time with population size N when mutational effects are weaker than 
the selection (regime I). The simulation data is obtained by averaging over 10 4 independent realisations 
for a = l/2(o, +), l(o, x), 2(A, *) for b = and 0.01/ respectively. The points shown by filled squares 
are the result of numerical integration of (|57p while the solid and dashed curves are the analytic 
formulae (|39p . The other parameters are / = 0.004, s = 0.2, U = 0.005 and A = 5. The right bottom 
inset shows the data for the same parameters as in the main figure for small populations and the left 
top inset shows the nonmutator fraction as a function of time for N = 128 and a = 1. 
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Figure 6: Variation of fixation time with population size N when mutational effects are stronger 
than the selection (regime II). The simulation data is obtained by averaging over 10 4 independent 
realisations for a = l/2(o, +), l(o, x),2(A,*) for 6 = and 0.01/ respectively. The filled symbols 
show the result of numerical integration of (|57p while the solid curves are the analytic formulae fj39fl . 
The other parameters are / = 0.0005, s = 0.001, U = 0.005 and A = 10. The inset shows the data for 
the same parameters as in the main figure for small populations. The simulation data for a = 2 does 
not match quantitatively with the diffusion theory possibly because k* = (U/s) 1 / ~ 6.7 is rather 
small. 



30 



10 4 




i 

o 


1 ^ 

o 

A 


lh- 


o 

© o 

o 


A 




10 3 


A 

A 
A 

i i 


i 


i 



50 100 150 200 

N 



Figure 7: Variation of fixation time with population size when rare beneficial mutations can also occur. 
The simulation data for regime I (/ = 0.004, s = 0.2, U = 0.005, A = 5, a = 1) are shown by triangles 
and regime II (/ = 0.0005, s = 0.001, U = 0.005, A = 10, a = 1) by circles. The blue points show the 
fixation time if a deleterious mutation occurs at a rate {7(1 — e) and a beneficial mutation at a rate Ue 
for e = 0.1. The red points show the simulation data in Figs. [5] and [6] for a = 1, b = for comparison. 
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